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ABSTRACT 

We present 880 jam images of the transition disk around the star [PZ99] J160421. 7-213028, a solar- 
mass star in the nearby Upper Scorpius association. With a resolution down to / /34, we resolve the 
inner hole in this disk, and via model fitting to the visibilities and spectral energy distribution we 
determine both the structure of the outer region and the presence of sparse dust within the cavity. 
The disk contains ~0.1 Mj up of mm-emitting grains, with an inner disk edge of about 70 AU. The 
inner cavity contains a small amount of dust with a depleted surface density in a region extending 
from about 20 - 70 AU. Taking into account prior observations indicating little to no stellar accretion, 
the lack of a binary companion, and the presence of dust near ~ 0.1 AU, we determine that the most 
likely mechanism for the formation of this inner hole is the presence of one or more giant planets. 
Subject headings: Stars: pre-main sequence; (Stars:) planetary systems: protoplanetary discs; sub- 
millimeter: stars 



1. INTRODUCTION 

It has been long suspected that planets form in disks 
of circumste llar material, with the hypothesis dating to 
|Kant| ( [T755| ). It was nearly two hundred fifty years, how- 
ever, before tools were developed to allow the direct ex- 
amination of this hypothesis. Large scale studies of ther- 
mal emission from small dust grains in circumstellar disks 
were enabled by t he Infrared Astronomical Satellite (e.g. 



Strom et al. 1989), while the development of millime- 
ter observatories allowed the direct measurement of dust 



mass ( Beckwith et al.|l990 ). The Spitzer space telescope 
allowed tor statistically complete censuses of nearby star 
forming regions (e.g. c2d, FEPS, Gould's Belt Survey), 
while ground based interferometers such as the Smithso- 
nian Sub-Millimeter Array (SMA) and Plateau de Bure 
Interferometer have allowed for imaging of t he dust and 
gas components of the brightest d isks (e.g. |Andrews fc| 



Williams|[2007 



Schaefer 



et al.||2009|) . 

also seen the 



Recent decades have also seen the discovery of hun- 
dreds of extrasolar planets. Ground-based radial ve- 
locity and transit searches have determined that ~10- 



20% of stars hos t Jupiter mass and larger planets (John 



son et al. 



20101). In the past year, the Kepler space 
telescope has revealed a similar fraction of st ars host- 
ing pl a nets d own to Neptune and Earth sizes ^Howard 



et al. (2011). Planet formation appears to be a com- 



mon process, though only in the past year has a planet 
been observed in a young disk, confirmin g the connec- 



tion betwe en disks and planet formation (Kraus & Ire- 
land||20lT ). The formation timescale remains uncertain, 
however. 

A likely sequence is that the disk sees the growth of 
dust grains, which aggregate to form planetesimals and, 
over several millions of years, planets. In the core accre- 
tion model, giant planets form in cases where the disk 
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still has sufficient gas mass for large, ~10 M e , plan- 
etesimals to experience gas accretion, and then runaway 
gro wth, a process t hought to occur in disks 3-10 Myr 



old 
of 



Pollack et al. 



19961). Whereas there are a handful 
Myr old massive disks, the medi an lifetime for in- 
frared excess from disks is ~3 Myr ( |Williams fc Cieza 
How ever , the recent millimete r surveys of |Lee 



2011) 



et al. (2011) and Mathews et al. (2012) have found only 
a handful of stars in the 1C 348 and Upper Scorpius re- 
gions hosting massive disks. Most disks are already well 
advanced in the process of evolution from gas and dust 
rich primordial disks to sparse dust debris disks. The 
paucity of massive disks in these regions compared with 
the relatively high frequency of exoplanets suggests that 
planet formation must be almost complete by the 3 to 5 
Myr age of these regions. 

In addition to placing constraints on the timescale of 
formation, observations of disks may provide clues to sys- 
tems where we have the best chance to directly observe 
planet formation in action. A gi ant planet is expected 
to tid ally open a gap in the disk (Artymowicz & Lubow 
1994), leading to inner disk dust depletion. Such deple- 



tion is observable as a reduction in near-IR emission or a 
hole in emission which could be imaged with high resolu- 
tion (sub) mm observations. Near to mid-infrared pho- 
tometry and spectroscopy with Spitzer has shown that 
less than ~10% of stars with disks show some degree of 
inner disk depletion yet still retain massive outer disks 



(e.g. Cieza et al. 2010). To date, about a dozen such 



transitional disks around single stars have had their in- 



ner d isk gaps verified by s ub-m illimeter imaging (jHughes 
|et al.|20Q7| |Brown et al.||2008| |Andrews et"a^201ip . 



Here, we present spatially resolved 880 /im contin- 
uum and CO J=3^2 images showing a transitional disk 
around a star in the nearb y Upper Scorpius asso ciation. 
At an estimated age of 5 ( jPreibiscfi et al.|[2002[ ), Upper 
Scorpius lies in the mid-range of the theoretical timescale 
for giant planet forma tion, and at a distance of 145 pc 
( de Zeeuw et al.p 999), it is one of the closest regions of 
recent star formation. Recently, a somewhat older ag e 
of ~12 Myr has been suggested by Pecaut et al. (2011), 
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which would place the association at the late-range of 
theoretical planet formation time-scales. 

The star [PZ99] J160421. 7-213028 (aka RXJ1604.3- 
2130A, henceforth J1604-2130) was first identified as 
a member of the Upper Scorpius OB association in 
the spectroscopic survey of X-ray selected sources of 
Preibisch fc Zinne cker| fl!999|). We adopt the stellar pa- 

sch e 

1 



rameters found by IPreibisch et al. (2002): spectral type 
K2, stellar mass of 1.0 ik/ , extinction of Ay 
log T eff of 3.658, and log L/L Q of -0.118. 

J1604 -2130 is included in the UCAC3 proper motion 
survey ( jZacharias et al.|2010[ ) . They report a J2000 posi- 
tion of a 20 oo =16:04:21.6541 and £ 2 ooo = -21:30:28.496, 
with errors on the J2000 position of 21 mas and 17 mas 
in right ascension and declination, respectively. They 
also report a proper motion of -14.5±1.9 mas/yr in right 
ascension, -17. Oil. 5 mas/yr in declination. 

Another association member, RXJ1604.3-2130B, lies 
16" to the southwest (proje cted separat ion 2350 A U) and 



is a potential companion ( Kraus & Hillenbrand 



Aperture masking inter ferometry by Kraus et al. 



2009). 

( pro ) 



ruled out the presence of a close binary companion of 
AK = 3.5 mag at a distance of 10-20 mas, and AK = 6.2 
mag at 40-80 mas (corresponding to — 0.06 M at a 
distance of -2 AU, and -0.01 M© at -9 AU), while 
direct imaging in that work ruled out a binary companion 
as faint as AK —5.5 at 240-320 mas (corresponding to — 
0.02 Mc?) at a distance of —4 AU). The complementary 



study of Ireland et al.| ( |2011[ ) has placed upper limits on 
companion masses of 0.07 M@ to 0.005 M at separations 
from 60 to >300 AU. 

High resolution optical spectroscopy shows evidence 
for very low levels of accretion, with the system exhibit- 
ing a P-Cygni profile in Ra emission (10% full width of 
252.7 km/s, equivalent width of -0.57 - -0.27 A, and he 



liocentric radial ve l ocity of -6.27 - -6 .90 km/s; Dahm 



& Carpenter|2009; 


Dahm et al. 


2011). Recent Spitzer 


1RAC, IRS, and M 


LIPS observations dCarpenter et al. 


2006 


Dahm & Carpenter 


2009 


Carpenter et al. 2009) 



and 8 //m, and a rising excess at 16 fim and longer wave- 
lengths. Dah m fc Carpent er 2009 noted this distinctive 



SED, as well as evidence of a factor of 4 variability in the 
near-infrared excess indicating possible changes in the 
dust environm ent in the inner disk on short timescales. 
|Dahm| ( |2010[ ) presented J, H, and K band spectroscopy 
of this target, with weak H and K band excess consistent 
with blackbody emission at 900 K, suggesting the pres- 
ence of dust at small orbit al radii (<0.1 AU). T he 1.2 



mm photometric survey of Mathews et al. (2012) shows 
that this disk is the most massive in Upper Scorpius. 

In ^2j we describe our (sub) mm observations and data 
reduction. We then show the continuum and CO line 
maps of the J1604-2130 transition disk, as well as anal- 
ysis of the visibility data sets and mass estimates, in ^3j 
Using a radiative transfer code, we model the dust com- 
ponent of the disk and estimate its parameters in Section 
[4] In ^5] we discuss the implications for disc evolution 
and planet formation, and show that the disk has likely 
formed one or more giant planets. We conclude with a 
summary in Section [6j 



2. OBSERVATIONS 

Interferometric observations of J1604-2130 were car- 
ried ou t in three trac ks with the Sub-Millimeter Array 
(SMA, |Ho et al.|2004| ) and in tw o tracks with the Plateau 
de Bure interferometer (PdBI, |Guilloteau et al. 1992). 
In Table [I] we give the night, array configuration, total 
continuum bandwidth, total time on-source, zenith sky 
opacity at 225 GHz, and median system temperature for 
each observation. 

SMA observations consist of upper and lower band- 
passes of 4 GHz bandwidth (2 GHz prior to 2009) for 
a total continuum bandwidth of 8 GHz (4 GHz prior to 
2009), with the centers of each bandpass separated by 5 
GHz and a local oscillator frequency of 340.4 GHz (880 
jam). Continuum regions were sampled at 3.25 MHz, 
and a 104 MHz wide region was set aside at the rest 
frequency of the CO J=3^2 line (345.79599 GHz) to 
carry out higher spectral resolution observations with a 
sampling of 0.203 MHz. Extended array observations in- 
cluded baselines from 44 to 226 m (50 - 257 kA), and very 
extended array observations included baselines from 120 
to 509 m (136 - 578 kA). Doppler tracking was carried 
out in the kinematic LSR rest frame. Our observations 
were carried out in generally good to excellent conditions, 
with system temperatures of —100 - 400 K, and zenith 
opacities at 225 GHz (T225) less than 0.1. 

For gain calibration, we interleaved alternate 4-5 
minute observations of the quasars QSO B1514-24 and 
QSO B 1730- 130 between 20-26 minute observations of 
J 1604-2 130. For bandpass calibration, we used the bright 
quasars 3C454.3 or 3C273, and Uranus, Callisto, or Ti- 
tan for flux amplitude calibration (as available). Obser- 
vations of flux and bandpass calibrators were carried out 
before or after J1604-2130 was available. 

Using standard routines in the facility IDL package 
MIR, we flagged and calibrated the data. We carried 
out baseline based phase calibration, finding rms phase 
errors of -10-20°. Based on variations in the measured 
fluxes of our gain calibrators, the flux calibration has an 
uncertainty of —15%. 

We combined wideband continuum channels from all 
observations and both sidebands in the imaging process. 
We used the facility reduction tool MIRIAD to carry 
out Fourier inversion, CLEAN deconvolution and image 
reconstruction using natural weighting. The synthesized 
beam of the combined observations is / /51x0 / /32, with an 
rms of 1.3 mJy/beam in the continuum. We separately 
carried out the same imaging process on the higher spec- 
tral resolution line data, binning in 0.25 km/s velocity 
channels centered on the CO J=3— >2 line. The beamsize 
is 0'/53 x / /34, and the median noise in the channels is 
125 mJy/beam. 

PdBI observations were carried out on two nights, Jan- 
uary 30, 2011 and March 2, 2011, in the A and B con- 
figurations, respectively. Continuum observations at 2.6 
mm used the WideX receiver with 3.6 GHz upper and 
lower sidebands for a total bandpass of 7.2 GHz and sam- 
ple spacing of 1.95 MHz. Line observations were carried 
out in a 20 MHz wide region at the rest frequency of 
the CO J=1^0 line (115.271 GHz), with a sampling of 
0.039 MHz. The A configuration observation included 
baselines from 136 to 760 m (52 - 292 kA), while the B 
configuration observation included baselines from 88 to 
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Table 1 

Observations 





Date 


Configuration 


Bandpass 




T 225, zenith 


median T sys 


2008, 


July 11 


SMA Extended, 345 GHz 


4 GHz 


2.0 hours 


0.05 


211 K 


2009, 


August 31 


SMA Extended, 345 GHz 


8 GHz 


2.7 hours 


0.09 


331 K 


2010, 


March 2 


SMA Very Extended, 345 GHz 


8 GHz 


5.8 hours 


0.03 


139 K 


2011, 


January 30 


PdBI A, 115 GHz 


7.2 GHz 


4.4 hours 


0.17 


466 K 


2011, 


March 2 


PdBI B, 115 GHz 


7.2 GHz 


4.9 hours 


0.12 


418 K 




2 1 0-1-2 1 0-1-2 1 0-1-2 

A RA (") A RA (") A RA (") 



(a) 880 fim continuum (b) Integrated CO 3-2 (c) CO peak velocity 

Figure 1. A) 880 /im continuum map, with contours indicating 2, 3, 6, 9, and 12 a intensity (lcr = 1.3 mJy/beam). B) Integrated CO 
J=3— >2 map from 4 to 5.5 km/s. Contours indicate 2, 3, 4, 5, 6, and 7a emission, in la steps (lcr = 75 mJy km/s /beam ). C) Peak CO 
J=3— >2 velocity map, limited to regions with 2a integrated line emission. In a) and b), the stellar position at 16:04:21.645 -21:30:28.66 
is indicated with a cross, and the solid and dashed lines in c) indicate the directions of the major and minor axes, respectively. The 
~ 0'/5 X 0'/3 beam is shown in the bottom left. 

452 m (34 - 174 kA). System temperatures ranged from 
400 to 800 K, and zenith r 22 5 ranged from 0.12 to 0.17 
(corresponding to median water vapor column from 2.3 
to 3.5 mm PWV). 

Gain calibration was carried out using 3-6 minute ob- 
servations of 1504-166, 1519-273, and 1657-261 (as avail- 
able) between 25 minute observations of J 1604-2 130, 
while flux and bandpass calibration were carried out with 
observations of MWC349 and 3C279. We used the facil- 
ity reduction tool GILDAS to carry out data calibration, 
Fourier inversion, and imaging. The combined contin- 
uum observations have a beam size of 2 / /01 x 0'/70, with 
an rms of 0.25 mJy/beam in continuum. From repeated 
observations of the calibrators, we estimate a flux cali- 
bration uncertainty of 20%. 

3. RESULTS 

We discuss first the appearance of the SMA maps, and 
then analyze the velocity map and continuum visibilities, 
along with continuum photometry from the literature, to 
determine the disk properties. In the PdBI maps, the 
disk is barely resolved in the east-west direction. As the 
SMA data is at much higher resolution, we use the PdBI 
continuum measurement as a point in the SED. From the 
continuum data, we measure the total flux at 880 fim and 
2.6 mm. From the integrated emission maps we measure 
the total flux in the CO J=3^2 and J=1^0 lines, and 
from the CO J=3— )>2 peak-velocity map we determine 
the orientation of the disk. We then use the continuum 
radially averaged visibilities to infer the central position 



Table 2 

Observational Properties 



Property Value 



F880/^m 164±6 mJy 

Fco J=3^2 5.2 ±0.1 Jy km/s 

F2600yLtm 5.1±0.5 mJy 

Fco J=i^o 0.48 ± 0.04 Jy km/s 

central position 16:04:21.645 -21:30:28.83 

inclination 6°±1.5° 

position angle -5° ±10° 

systemic velocity (LSRK) 4.7 km/s ±0.1 km/s 



and the inner radius of the disk. These parameters are 
then refined via modeling, in Section HI We list the prop- 
erties derived in this section in Table[2j 

3.1. SMA images 

We show the 880 am continuum emission map for 
J1604-2130 in Figure [l] centered on the disk central po- 
sition. The disk is being viewed from a nearly face-on 
perspective, and has a large inner hole. The image has a 
peak flux of 18.2 mJy/beam. To measure the total emis- 
sion, we integrate all flux above the 2a contour, giving a 
total flux of 165±6 mJy. 

CO J=3^2 line emission was found at 2a and greater 
significance in channels from 4.0 and 5.5 km/s (LSRK). 
We sum these 0.25 km/s wide channels to produce the in- 
tegrated intensity map seen in Figure [T] with an rms of 75 
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Figure 2. Real flux components of the continuum and CO J=3— >2 visibil ities, la er rors are shown as vertical bars, and model fits (i.e. 
integrated zeroth order Bessel functions of the first kind, following Hughes et al. 2007) shown as dashed red lines. The range of 3a values 
for the imaginary components are shown as the gray region. The position ot tne hrst null strongly constrains the radius of the inner edge 
of the circumstellar disk (50 - 90 AU, for varying radial emission functions). 



mJy km/s / beam. The integrated intensity map shows 
a drop in emission to the southwest and in the center 
of the disk; however, the emission does not clearly drop 
below 3<j in either location. The outer radius of the CO 
emission roughly corresponds to that of the continuum 
emission. Integrating emission above the 2a contours of 
integrated line emission, we find a total CO J=3^2 line 
flux of 5.2 ± 0.1 Jy km/s. Following the same procedure, 
we find a line flux of 0.48 ± 0.04 Jy km/s for the CO 
J=1^0 data. 

We generated the peak velocity map by fitting a gaus- 
sian emission profile at each spatial position in the disk, 
including emission from 0.75 to 8.5 km/s in order to 
achieve good baseline fits and ensure we captured the 
full profile at each location. The median uncertainty in 
the peak velocity at each point is 0.07 km/s. This map 
shows a clear east- west velocity gradient which we use to 
measure the position angle and inclination of the disk. 
We do this by fitting a Keplerian velocity profile (Fig- 
ure [I]). Velocities along the observed minor axis of the 
disk should be constant, reflecting the systemic velocity. 
Velocities along the major axis will be proportional to 
r _1 / 2 sim, where r is the distance from the star and i is 
the inclination from our line of sight (with i = indicat- 
ing a face-on disk). 

To deduce the observational geometry, we varied the 
central position, position angle, and inclination, calcu- 
lated line-of-sight velocities along the major and minor 
axes, and carried out a x 2 fit to the velocity map. We 
mask the central 0'/53 (the beam major axis) of the ve- 
locity map, since beam dilution will mask the large range 
in velocities expected at small radii and lead to poor fits. 
We find the disk is centered at a^ooo = 16h04m21. 6416s 
and £2000 = — 21°30'28 / /798, with a conservative uncer- 
tainty of ~ 0'/l in both right ascension and declination. 
The predicted position from the UCAC3 position and 
proper motion is 16h04m21. 6394s, -21°30'28'/661. Tak- 
ing into account the J2000 position uncertainty, proper 
motion uncertainty, and the position uncertainty of our 



observations, the offset between the predicted stellar po- 
sition and the disk center is Aa = 0V031 ± / /104 and 
A5 = / /137±0 / /102, respectively. The disk is inclined at 
6 ± 1.5° from face on, at a position angle of 5° ±10° west 
of north, and with a central velocity of 4.7 ± 0.1 km/s. 

We use these parameters for modeling the dust disk, 
but defer further quantitative discussion of the line data 
to a later paper. 

3.2. Radially averaged visibilities 

Based on the disk inclination and position angle as de- 
termined from the peak velocity map, we deproject the 
interferometer visibilities, scaling uv distances along the 
minor axis by the cosine of the disk inclination. We plot 
the radially averaged real and imaginary fluxes versus 
the radial uv distance. To provide an adequate signal- 
to-noise ratio, we bin the real and imaginary flux in 25 
kA bins from 12. 5-287. 5kA, and 50kA bins at larger base- 
lines. 

We plot the binned components of the continuum vis- 
ibilities in Figure |2j along with the similarly binned CO 
J=3^2 visibilities, including measurements from 4 to 5.5 
km/s as in the integrated emission map. The imaginary 
component is less than 3<r in all bins, making it consis- 
tent with zero as expected for a centered, azimuthally 
symmetric disk. We find that a greater fraction of the 
power in the CO line is found at small baselines com- 
pared with the continuum emission, suggesting more CO 
emission is to be found in large scale structures. 

The radially averaged profiles strongly constrain the 
size of the inner hole in gas and dust, and provides further 
insights to the large scale structure of the disk. Strong 
contrast between the cavity and outer disk leads to nulls 
in the Fourier transform of the emission profile. To esti- 
mate the inner and outer edges of sub-millimeter emis- 
sion, we fit a simple disk model to the ra dially binned 
real v isibilities based on Equation A2 of Hug hes et al. 
(pQ07|>: 
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Vis(R uv ) = A [ e 1 - p J (27r9R uv ) d0., 



(1) 



which represents the Fourier transform of a power-law 
brightness disk with an inner hole (i.e. 1(6) oc 0~ p ). A is 
a free parameter that sets the intensity scale, 6\ and 62 
are the inner and outer disk edges respectively, and p is 
the exponent of a power law describing the variation in 
emission intensity with radius. Jo is the Bessel function 
of the first kind of order 0. We show our best-fit analytic 
models in Figure [2] The continuum model has parameter 
values of A = 0.16 ± 0.02 Jy, t = 0743 ± 0701, 2 = 
0795 ± 0703, and p = 2.03 ± 0.01. The CO line model 
has parameter values of A = 5.0 ± 0.2 Jy km/s, 0\ = 
0732 ± 0701, 6> 2 = 1784 ± 0724, and p = 1.98 ± 0.01. The 
similar values for p (~ 2) suggest the brightness functions 
of dust and gas have similar variation with radius. 

The inner and outer angular limits correspond to radii 
from 60 - 140 AU for continuum emission, and 45 - 270 
AU for CO emission. These suggest the dust cavity may 
be larger than that of the gas, however, given the com- 
plex interaction of surface density and temperature struc- 
ture in determining the surface brightness function, such 
an assessment cannot be made without more detailed 
modeling. We also note that the continuum model does 
not fit the long baseline visibilities well, suggesting the 
presence of some millimeter continuum emission at small 
spatial scales. We examine the structure of the dust disk 
with detailed radiative transfer modeling, including the 
spectral energy distribution, but we defer further inter- 
pretation of the CO line data to a later paper. 

3.3. Millimeter- emitting dust mass 

Under the assumption of optically thin dust emission, 
we can estimate the dust mass of the disk. The dust mass 
(Mdust) wi ll be directly propo rtional to the sub-mm flux 
(F„), 



as m 



Hildebrand (1983): 



M du 



d 2 F v 
n v B v {T c ) 



(2) 



We assume J1604-2130 lies at the mean distance of 
Upper Sco derived from the Hipparcos data (145 pc, de 
Zeeuw et al. 1999). We also use a characteristic dust 
temperature (T c ) of 20 K as shown to be appropriate 
for Taurus disks ( Andrews fc Williams|2005] ), and a dust 
mass opa city k v = 10 (z//1000 GHz)^ cm 2 /g ( jBeckwith] 
eTaI1[l990l ). 

We estimate the value of /3 by making a power law fit 
to the mm- fluxes. We plot the 880 /im and 2.6 mm fluxes 
in Fig ure [51 along with the 1.2 mm flux from |Mathews| 
et al.| ( [2012| , 67.5 ± 1.4 mJy. A power law fit, including 
calibration uncertainties, shows that for f v oc u~ amm 5 
the mm-slope a mrn is equal to 3.17 ±0.21. Assuming 
an optically thin medium and emission in the Rayleigh- 
Jeans limit, where a mm = /?+2, we then calculate a value 
of (3 = 1.17. This is approximately 3a less than the value 
of (3 for small ISM grains (/? = 2), and is greater than 
that expected for large blackbody grains (/? = 0). It is 
close to the value of f3 ~ 1 expected in the case where 
grains follow a power law size distribution N(a) oc a -3,5 , 
with a maxi mum grain siz e at least 3 times the observed 
wavelength QDraine| |2006 ) . Together, these suggest that 
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Figure 3. We show the millimeter SED as black dots, along with 
power law fit and comparisons to Q; mm = 2 and 4, corresponding to 
opacity power laws of {3 = and {3 = 2 in the case of optically thin 
emission. Errors include both formal and calibration uncertainty. 

this disk may have experienced significant grain growth, 
up to at least centimeter-sizes. The inferred dust mass 
is 0.1 M Jup . 

The assumption of optically thin emission is justified 
because, for the f3 = 1.17 opacity function, dust will be- 
come optically thick to 880 /im emission at a surface den- 
sity of 0.35 g/cm 2 , and at a surface density of 1.3 g/cm 2 
for 2.6 mm emission. Surface densities this high are not 
observed at large radii in even the youngest primordial 
disks. For our estimated disk mass, this high density 
would require a 2-5 AU wide ring at a radius of 90 AU. 
However, such a narrow ring is ruled out by both the 
SMA images and the simple model fits describ ed a bove. 
In addition, we show in our later modeling (j ]4.2| ) that 
the peak surface density in the disk is a factor often too 
low for the dust to be optically thick to 880 jam emission. 

4. MODELING 

The temperature profile of the disk is reflected in the 
spectral energy distribution (SED), while the density 
profile is reflected in the resolved sub-mm image. By 
simultaneously fitting the SED and SMA data, we de- 
rive the disk properties. Because of the sparse sampling 
of the Fourier plane, we fit the interferometric visibilities 
rather than the image itself. For the stellar properties, 
we adopt a stellar temperature of 4550 K, luminosity of 
0.76 L , logg of 4.0, and extinction of 1 mag in V band, 
based on the references in ^T] For the disk orientation, 
we use the d isk position angle and inclination derived in 
Section 13.11 

4.1. Procedure 

To model the disk, we use the 3-dimensional r adiative 
transfer code MCFOST flPinte et al.||2006| |2009[ ) to pro- 
duce a model SED and 880 /im image. We then simulate 
observation of that image to generate visibilities. We fit 
the SED and radially binned components of the 880 /im 
continuum visibilities. 

We constructed a spectral energy distribution for 
J1604-2130, combining our 880 /im flux with values in 
the literature. For the B, R, and I bands, we used 
photometry from Preibisch & Zinnecker (1999), convert- 
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ing magnitudes to flux es using the absolute flux conver- 
sions of Bessell] (|1979|. We used near-infrared J, H, and 



band ph otometry from the 2MASS survey (Sknitskie 



et al.l 2006]) and the absolute flux conversions ol |Cohen 
et al.| (|2UU3|). The near-infrared 4.5, 8, and 16 /im pho 



tometry ol |Carpenter et aL (|2QQ6) and mid- infrar ed 24 
and 70 /im photometry ol |Carpenter et al. (2009) were 
published as fluxes; these values are used here directly. 
We also include the spatially integrated 880 /im flux from 
this work, and the single-dish 1 .2 mm continuum mea- 
surement from |Mat hews et al.| (2012). We dereddened 
the optical and near-infrared fl uxes by Ay = 1 m ag us- 
ing the R=3.1 reddening law of Fitzpatrick (1999). 

The SMA image indicates a disk with a large inner 
hole. In addition, the SED shows only photospheric emis- 
sion out to 8 /im, requiring an absence of dust out to sev- 
eral AU. However, the excess emission at 16 and 24 /xm 
requires the presence of a small amount of dust within 
the inner hole seen in the SMA image. A model that has 
an inner hole entirely depleted of dust therefore does not 
fit the data, and we require a thr ee zone disk model. 

Following | Andrews et al.| ( |2011[ ), we explore a modified 
version of tne parametric accretion disk density model 
with a dust surface density of 



D \ -7 

^ acc (R) = S C ( — J exp 



RcJ 



(3) 



where £c is the surface density at the characteristic scal- 
ing radius, Rc, and 7 is treated as a free parameter de- 
scribing the surface density of the disk (this corresponds 
to t he viscosity power law index in accretion disk theory, 



e.g. 



Hartmann 



We split the 



et al.||1998[ ). 

disk into thi 



ree regions. In the innermost 



part of the disk, within a gap of radius R gap , we treat 
the dust surface density as zero to account for the lack of 
excess at wavelengths up to 8 jam. Based on the hole seen 
in 880 /im emission, yet accounting for the excess seen at 
16 and 24 /im, we deplete the inner region of the disk by a 
factor S cav from R gap to R ca v We refer to this depleted 
region as the disk cavity, and to the remainder as the 
disk. For the outermost radius included in our models, 
we set Rout to 300 AU, but in practice, the emission 
at these radii is too low for the sensitivity of our data. 
Therefore, our adopted surface density model goes as: 



^dust 



(R) 



$cav ^acc 

s, 



0, 
OR), 
OR), 



if R < Rg a p 

< R < R r 



if i£ 

if Rr.av < R < R n 



gap 



(4) 



The vertical dust distribution is modeled as a Gaussian 
with scale height h = hcRc (R/Rc)^ , where hcRc is 
the scale height at Rq and i/j describes the power-law 
disk flaring. After initial modeling indicated our data 
do not constrain -0, we adopted a value of 1.2 for this 
parameter, based on the median flaring value found in 



the multi-transition disk study of | Andrews et al. (2011). 
For the dust properties, we adopt the porous grain min- 



eralogy of |Mathis fc Whiffen (1989), fixing the exponent 
of the power law grain size distribution at —3.5. We 
adopt a minimum grain size of 0.005 /im and a maxi- 
mum grain size of 7.8 mm (i. e. 3 times th e longest ob- 
served wavelength, following |Draine 2006). The grain 



composition and grain size distribution determine the 



wavelength-dependent opacity function. We also no te 
that the choice to use the grain properties of Mathis 



Whiffen is essentially arbitrary. Other grain models lead 
to equally good models and the current data on J 1604 
does not discriminate among the several possibilities. 

In order to find a best-fit model for the J1604-2130 
disk, we carried out a search of the 7 free parameters 
identified above (£<?, Rc, 1, Rgap, Rcav, 5 cav , he), 
adopting the central position, inclination, and po sition 
angle determined from the peak- velocity map (j |3.1| ). 
Generating even a sparse grid of 7 parameters would 
require a prohibitively long time, therefore, we use the 
Levenberg-Marquardt x 2 minimizatio n algorithm as im- 
plemented in the IDL routine MPFIT (Ma rkwardt|2009 ) 
to carry out an iterative least squares lit to the observed 
SED and real and imaginary components of the radially 
averaged visibilities. Each evaluation of the parameter 
search carries out the following steps: 



1. For a given set of parameters, we generate the 
corresponding model at the observed central po- 
sition, position angle, and inclination. MCFOST 
produces both an SED for the model and an 880 
/xm image. 

2. We then simulate an observation of the model with 
the uv-coverage and integration time of our SMA 
observations. We use the Fourier transform func- 
tionality of the CAS A packagqjto generate a sim- 
ulated uv dataset for the modelimage, from which 
we generate deprojected, recentered, and binned 
radially averaged visibilities. 

3. The fitting routine calculates the local numerical 
gradient of the \ 2 function of the SED and the ra- 
dially averaged visibilities. As the model disk is az- 
imuthally symmetric, the imaginary components of 
the radially averaged model visibilities are always 
zero, with slight variations due to the stochastic 
nature of the image simulation. 

4. The numerical gradient of the x 2 values are then 
used in determining the next point to sample in 
the space of free parameters. The search algorithm 
follows the x 2 gradient towards the minimum. 

In order to avoid local minima and more fully map 
the parameter space, we carried out the search algo- 
rithm 40 times, randomly varying the starting values of 
the parameters within the ranges listed in Table [3j The 
Levenberg-Marquardt algorithm naturally leads to low 
density sampling in regions of the parameter space with 
high x 2 values, with the sampling density rising as the 
search moves into regions of lower x 2 • In addition to pro- 
viding a sampling of the parameter space, each run, fc, 
of the search algorithm results in a set of best-fit values 
for the i parameters (/i^ = /ii,fc, /i2,fc 5 •-■,/ J <i,k), parameter 
uncertainty estimates (dk = <Ji,k,&2,k, •••?c r i,/c) 7 and a co- 
variance matrix (Ak) which allows for estimation of the 
parameter dependencies. 



3 http://casa.nrao.edu/ 
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Table 3 

Dust parameters for USco J 1604 



Parameter 



Mean Standard Search 

value deviation Lower Upper 



log(S c ) [g/cm' z ] 


-1.23 


0.17 


-3 


1 


Rc [AU] 


106 


7 


90 


150 


7 


-0.40 


0.43 


-2.0 


2.0 


Rgap [AU] 


20 


3 


15 


49 


Rcav [AU] 


72 


3 


50 


100 


\og(5 C av) 


-0.86 


0.35 


-3 





he 


0.04 


0.01 


0.010 


0.100 


M dust (M Jup ) 


0.13 









4.2. Fit Results 

We carried out 40 parameter searches starting from 
very different initial values, each of which consisted of 
about 100 disk models. The model with the lowest 
X 2 value out of the total of 3967 separate disk models 
has x 2 = 75.7 and parameters of log Sc a v — —0.89, 
R c = 108AU, 7 = -0.69, logS c = -1.23 g/cm 2 , 
h c = 0.04, R gap = 19 AU, and R cav = 70 AU. 35 of 
the 40 searches converged to a very similar set of pa- 
rameters, with x 2 ranging from 76 to 85. This gives us 
confidence that we are reaching the global x 2 minimum 
of the seven dimensional parameter space. The remain- 
ing 5 of the 40 searches converged on local minima with 
X 2 greater than 100. 

The separate searches show strong convergence, de- 
spite starting in different regions of the parameter space. 
These models tend to have similar x 2 values, suggest- 
ing that a mean of the parameters found by the several 
searches will provide a reasonable estimate of the disk 
parameters. We make use of the full set of ~ 4000 disk 
models to map out parameter space so we can assess the 
reliability of our parameter estimates. 

There are, however, difficulties in estimation of the pa- 
rameter uncertainties. Formally, our models have 40 de- 
grees of freedom (i.e. 14 SED points, 17 visibility bins 
with both real and imaginary components, and 7 free pa- 
rameters), for a minimum reduced x 2 value of ~ 2. In ad- 
dition, some data points do not contribute to the fitting 
process (e.g. optical photometry, where the disk cannot 
have emission, and imaginary visibility points which will 
be near zero for all azimuthally symmetric models), and 
others will be highly correlated (e.g. photometry near 
the SED peak and visibilities near the first null, both of 
which are strongly affected by the inner disk radius) and 
thus will not contribute a full degree of freedom. The 
combination of these two issues - possible model mis- 
specification and uncertainty in the number of degrees of 
freedom - suggests that the uncertainty estimates pro- 
duced by MPFIT for each search and the conversion of 
X 2 values to probabilities are inaccurate. 

This means that traditional approaches to error esti- 
mation (e.g. Ax 2 contours) are not applicable. However, 
the relative comparison of models via differences in x 2 
are still valid, i.e. a lower x 2 value still indicates a bet- 
ter model. We cannot say, however, to what degree the 
model is better. With these limitations in mind, we re- 
port the standard deviation of best-fit parameter values, 
but note that they do not necessarily represent la uncer- 
tainties, and only serve to give a qualitative indication 



of how well each parameter is constrained in relation to 
the range of values included in the search. We report the 
mean parameter values and their standard deviations, 
along with the range of values explored for each param- 
eter, in Table [3j 

For a further qualitative examination, we show the 2- 
dimensional A% 2 distributions for pairs of parameters 
(Fig. [4|. At each point within the displayed space, we 
calculate the parameter-distance and x 2 ~ we ighted mean 
X 2 - In each figure, we show best-fit parameters from 
our searches. The inner disk dust depletion factor, 5 cavi 
strongly correlates with both Rc and 7, the character- 
istic radius and surface density power law term. The 
other parameters (Ecs ^cs Rcav, an d Rgap) are tightly 
constrained with only weak correlations with other pa- 
rameters. 

We show the comparison of our adopted disk model 
with the data in Figures [5] (SED and binned visibilities) 
and [6] (continuum image, model, and residual). We cal- 
culate the residual by subtracting the model visibilities 
from the observed visibilities, and then carrying out the 
same imaging procedure as described in ^2] We can see 
that the model underpredicts the flux at the 2a level at 
the location of the peak of emission to the southeast, sug- 
gesting a possible asymmetry in the continuum emission. 

In Figure [7j we show the dust surface density for our 
adopted disk parameters. The total inferred dust mass 
is ~ 0.1 Mj up , in cl ose a greement with the simple calcu- 
lation presented in £3.3 The inner cavity in this model 
has a dust mass of ~ 0.007 Mj up . 





(b) 





Figure 4. 2-dimensional projections of the 7-dimensional x 2 dis- 
tribution ( 5cav appears twice, to highlight the strong covariance 
with Rc and 7). Gray circles show the results of the searches for 
the bestfit parameters and the black circle shows the parameters 
of the minimum-x 2 model. The apparent blocky structures in the 
high Ax 2 regions are artifacts of sparse sampling and high x 2 gra- 
dients. The projection of the low Ax 2 region to the left in Figure 
b is an artifact of the distance contribution term, as the minimum 
X 2 model contributes disproportionately to this nearby high x 2 
region, which has a low sampling density. 
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Figure 5. Model fits to SED and continuum visibilities. A) Spectral energy distribution for J1604-2130. We show our 880 /im and 2.6 
mm fluxes as large black circles, and dereddened optical, 2MASS, Spitzer, and 1.2 mm IRAM 30m observations as solid black circles. We 
assume calibration errors of 20% for optical photometry, 10% for infrared photometry, 15% for 880/im and 1.2mm photometry, and 20% for 
2.6mm phot ometry. Original photometry at short wavelengths is shown with black squares. The closest match stellar model from [Kurucz 
(1979 1993) is shown as a gray line. B) Real flux components of the 880 /im continuum as a function of deprojected uv distance, la errors 
are shown as vertical bars. 



5. DISCUSSION 

5.1. Disk mass 

Our continuum observations and modeling allow us to 
determine the dust mass in the disk of J1604-2130, and 
the high disk luminosity rules out the possibility of this 
being a debris disk. In addition, the detections of the CO 
J=3^2 and J=1^0 lines indicate the presence of gas in 
the disk. However, CO attains high optical depth at low 
column densities and provides only a lower limit to the 
gas mass. A more accurate determination requires multi 



with IRS spectroscopy which shows no evidence for th e 



10 or 20 /im silicate features (Dahm & Carpenter 2009), 
suggesting a deficiency of micron sized grains at ^1 AU 
scales. However, there are indications that the inner re- 
gion is not entirel y devo i d of m aterial. Near-infrared 
spectroscopy from Dahm (2010) reveals excess emission 
corresponding to black body em ission at a temperature 
of -900 K. From|Wyattl([2008l), 



rdust(AU) 



line s pectroscopy and detailed modeling (e.g .|Kamp et al. 
|2011| and is the subject of a future paper. To provide a 
rough estimate for the disk gas mass, we selec t from the 
DENT grid of disk models flWoitke et al.|2QlQ| ) all young 
(age <10 Myr), solar mass stars with lace on 0.1 Mj up 
dust disks and large inner holes. Among these, the detec- 
tion of CO J=3^2 at a level of - 10~ 19 W/m 2 suggests 
a gas mass of 3 Mj up with a factor of 10 uncertainty. 

5.2. Inner disk 

Our adopted model indicates the inner 20-72 AU of the 
J1604-2130 disk has 0.007 Mj up of dust in grains up to 
~1 cm in size. This represents a decrease in surface den- 
sity of a factor of ~10 relative to the dust rich outer disk 
shown in our SMA image. Examination of the parameter 
space indicates that all feasible models achieve a similar 
low inner disk dust density, either through this abrupt 
decrease in surface density or a combination of greater 
value of the characteristic radius and a much steeper de- 
cline in surface density (Rc as high as ~ 115 AU and 7 
as low as ~ —1.5, respectively). However, even this ex- 
treme requires a factor of 2.5 decrease in surface density 
at R C av The dust surface density decreases towards the 
center of the disk, and the transition cannot be smooth. 

The additional truncation of dust near 20 AU is indi- 
cated by the rising SED at 16 and 24 /im and is consistent 



278.3 



L 



0.5 



(5) 



this corresponds to blackbody dust orbitin g at approxi- 
mately 0.08 AU. |Dahm fc Carpenter| Q2009| ) noted a fac- 
tor of 4 variation between tiuxes measured at 4.5 and 8 
/im with Spitzer IRAC, and measurements at those wave- 
lengths a few weeks later with IRS spectroscopy; these 
wavelengths roughly correspond to blackbody tempera- 
tures between ~1100 and 600 K, which could in turn 
arise from dust at radii from ~0.05 to 0.17 AU. 

The CO imaging presented here does not rule out the 
presence of gas within the inner disk region. Indeed, the 
inverse P -Cygni profile and 10% full- width of 252 km/s 



ter 



Dahm et al. (2011) suggest t he system could be 

Dahm fc Carpen^] 



acc r eting, albe it at a very l ow ra te 



(2009) and Dahm et al. (2011) report Ra equivalent 
widths of -0.57A and -0.27A, respectively; this is too low 
for the system to be clearly identifi ed as accreting by the 



spectral type dependent criteria of |White fc Basri ( 2003 ) 
(3 A for K0 - K5 stars). The syst em similarly shows no 
emission in the K-band Br7 line ( |Dahm|[2010[ ), another 
accretion indicator. Usin g the relation s between accre- 
tion and Ha luminosity of 



Dahm 



upper limit on the accretion rate of 10 M /yr. 



(2008), we estimate an 
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Figure 6. A) 880 /im continuum map, with contours indicating 2, 3, 6, 9, and 12a intensity. B) 880 /im model image, with contours 
indicating 2, 3, 6, 9, and 12a intensity. C) 880 /im residual map, with solid contours indicating 2a residuals, and the dashed line indicating 
the 2a contour of the observed continuum map. In all images, the stellar position at 16:04:21.645 -21:30:28.66 is indicated with a cross, 
and la = 1.3 mJy/beam. 



5.3. Hole formation mechanism 

Several mechanisms could lead to the formation of 
an inner hole in a circu mstellar disk, summarized in 
Williams fc Cieza] ( |2Q11[ ). Recent work has provided a 
step- by-step classification scheme to identify the domi- 
nant mechanism leading to the formation of a transition 



disk. Cieza et ah] pOlOj) studied a large sample of young 
transition disks and developed a systematic classification 
system. Below, we discuss these potential mechanisms 
and address their feasibility in explaining the formation 
of the J1604-2130 transition disk. 

Grain growth: The mm slope of a mm =3.17 between 
880 /im and 2.6 mm indicates that substantial grain 
growth has occurred in the outer disk. In the warm inner 
regions of the disk, however, grain growth and dust set- 
tling are expected to lead to disks wit h a smooth down - 
ward slope to their infrared excess (e.g. |Cieza et al.|2010[ ), 
in contrast to the essentially non-existent excess shown 
by J1604-2130 at wavelengths up to 16 /im, where the 
SED abruptly rises. 

In addition, accretion is expected to continue in disks 
where an opacity hole has for med due to grain g rowth . 
From the empirical model of |Hartmann et al.| ([1998), 
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Figure 7. We show the surface density function for our adopted 
disk model. The surface density within the cavity (20 - 72 AU) is 
reduced by a factor ~ 10 relative to the disk. 



we would expect a 3 Mj up disk at 5 Myr to have an 
accretion rate ~ 3 x lO _lo M /yr, though this relation 
has large scatter for individual objects. This is more 
than an order of magnitude higher than the estimated 
upper limit of 10- n M o /yr for J1604-2130. While the 
accretion rate could be consistent with grain growth, the 
sharp changes in dust surface density suggest that while 
growth has occurred, it is not likely the primary source 
of the disk cavity seen in J1604-2130. 

Photoevaporation: As a disk depletes its gas con- 
tent, the accretion rate is expected to drop as well. When 
the accretion rate is comparable to the photoevapora- 
tion rate (typically ~ 10 -10 MQ/yr), photoevaporation 
is expected to become the dominant evolution mecha- 
nism for the disk, halt accretion, and clear the disk from 
the inside-out. This is expected to take place when disks 
have a total mass of 0.05-0.5 Mj up (I Alexander et al. 
2006). While this is approximately an order ot mag- 
nitude less than the mass of the J1604-2130 disk, the 
large uncertainty in the disk mass prevents us from rul- 
ing out photoevaporation as the driver of hole formation 
in this disk. However, it is unclear that photoevaporation 
can produce inner regions of depleted but non-zero dust 
mass. Such a region is not found in current models of 
disk structure under the influence of photoevaporation. 
Moreover, while the observed Ha equivalent width is too 
low for the system to be formally described as accreting, 
the inverse P-Cygni profile seen in Ha emission (Dahm 
et al.||20lT| ) suggests that accretion is continuing, albeit 
at too low a level to be precisely determined. 

Binarity: Disks can form holes due to tidal trunca- 
tion by a b inary companion (e.g. LkCa 4, Ireland & 
Kraus 2008| )). An equal mass companion opens a disk 
gap with a radius approxi mately twice the semi-majo r 
axis of the binary system ( |Artymowicz fc Lubowp 994 ) . 
However, at 36 AU the presence ot a binary companion 
to J 1604- 2 13 has been ruled out to an upper limit of 
~20 Mj up ( Kraus et al.||2008 ). Combining the imaging 
and apertur e-mas king mterterometry surveys of [Kraus 
et al. ( |2008| ) and |Ir eland et al.] (2011), a companion is 
ruled out at all radii within the disk down to ~13-20 
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Mj up , except for 2 AU and less (< 60 Mj up ), and ~60- 
75 AU (< 70 Mj up ). However, the limit on companion 
mass in the ~60-75 AU range is likely much lower than 
reported, as J1604-2130 would not have exhibited clo- 



sure p hases consistent with a point source in Kraus et al 



(|2008| ), and would have triggered further investigation 

(A. Kraus, personal co mmunication) . 

Planet formation: Cieza et al. (2010) use the lack 
of a near-infrared excess, followed by a steeply rising ex- 
cess, as an indicator for dynamical clearing. The SED 
of J1604-2130, which rises steeply from 16 to 24 /im, 
matches this criterion. However, in their study all sys- 
tems identified as potentially hosting a planet were still 
actively accreting. Given that a Jupiter mass planet can 
inter cept up to 90% of the gas accreting from the outer 
disk flLubow fc D'Angelo| [2006), a more massive planet, 
or several planets, could possibly be reducing the accre- 
tion rate onto the star. 

A single planet will directly clear only a narrow gap in 
a disk. However, after gap formation, inner material will 
accrete onto the star and leave a large gap. A typical 
accretion disk model with an evolution time of 1 Myr at 
30 AU experiences sig nificant evolution on a timescale 
of ~ 5 Myr at 70 AU (|Hartmann et aL|[l998| |Armitage 



2011), suggesting that a single giant planet could have 



opened the hole in the J1604-2130 disk. Disk cl earing ef- 
ficiency increases with a dditional giant planets (Do dson-| 
Robinson fc Sal^[2QlTl ). 

We have ruled out binarity as the mechanism for the 
formation of the gap in the J1604-2130 disk, and con- 
sider grain growth an unlikely candidate due to the two- 
region dust surface density function. The complicated 
dust distribution, in conjunction with marginal evidence 
for accretion, makes photoevaporation an unlikely cause 
for the inner hole. In our further discussion, we explore 
possible scenarios for a nascent planetary system. 

We can use the upper limit on accretion to estimate a 
limit to the size of a single giant planet in this system. 
Models explore the possibility of a planet truncating the 
disk, while allowing sparse m aterial to move past (e.g. 
Alexander fc Armitage||2007| . In such a scenario, the 
putative planet would be located close to the disk trun- 
cation radius; for illustration, we examine the case of a 
planet located just inside R ca v> at an orbital radius of (65 
AU) . The final truncation of the disk appears to occur at 
a smaller orbital radius, R o an, an issue we discuss below. 

[Alexander fc Armitage| ( |2007[ ) derived an expression 
relating the stellar accretion rate of a protoplanetary sys- 
tem to the disk mass and planetary properties. We solve 
their Equation 29 to estimate the mass of a putative 
planet, assuming the stellar accretion rate is approxi- 
mately 1/4 the planetary accretion rate: 
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with disk mass M^, accretion rate M acc , accretion pa- 
rameter a acc , orbital frequency ^, scale height at the 
planetary orbit planet orbital radius i?, and disk in- 
ner and outer radii R in and R ou t- From our disk model, 
we adopt Ri n and R out of 72 and 200 AU, respectively. 
The dimensionless accretion parameter a acc can reason- 
ably have values ranging from 0.001 to 0.01. 
From the Ra equivalent width of -0.57 A, we estimated 



an upper limit to stellar accretion of lO _11 M0/yr. At an 
orbital radius of 65 AU around a 1M star, the orbital 
frequency will be about 0.002 yr~ x , and assuming the 
100:1 gas-to-dust mass ratio of the ISM, the disk has a 
mass Md = 10 -2 M@. Our adopted disk model gives an 
estimated scale height, H, at 65 AU of less than 3 AU. 

Given these assumptions, for an accretion constant 
®acc = 0.01, a planet at 65 AU with a mass of 6.4 Mj up 
could allow accretion onto the central star at the upper 
limit rate of lO _11 M0/yr. With an increasing flaring 
constant, ip, and hence smaller scale height, the necessary 
planet mass decreases, reaching 6 Mj up for ip = 1.30. A 
lower value of ip would lead to a correspondingly higher 
planet mass, e.g. about 8 Mj up for i/j = 1.05. 

By comparison, a lower accretion rate could suggest 
a more massive planet; however, the observation of Ha 
emission in a P-Cygni profile suggests that while stellar 
accretion may be low, it is unlikely to be significantly 
lower than this upper limit. This in turn suggests the 
mass of such a planet could not be more than a factor of 
a few larger than 8Mj up , else accretion would be so low 
as to eliminate even the hint of Ha emission seen in this 
system. In addition, the planet mass scales with a acc . 
For a acc = 0.001, the planet mass could be an order of 
magnitude lower. 

We have so far ignored the trunca tion of the sparse 
inner disk at 20 AU. At this radius, Kraus et al. (2008) 
place an upper limit of 13 Mj up on the mass ol a com- 
panion. However, our disk model suggests a low dust 
density, which could be diverted by a low mass planet. 
This planet would be in addition to the mechanism caus- 
ing the decreased dust surface density within 72 AU. 

The complex dust structure seen in the J1604-2130 disk 
- variable dust at a fraction of an AU, a fully depleted 
region up to 20 AU, and a partially depleted region from 
20 to 72 AU - m ay be best explain ed by the presence of 
multiple planets. Zhu et al. ( 2011[ ) present a scenario in 
which four planets lead to a similar configuration, though 
they must invoke an additional du st filtering mechanism. 
Dodson-Robinson fc Salyk] ( |2011[ ) explore similar multi- 
planet scenarios, where low infrared opacities are driven 
by geometric filling effects as accretion continues through 
dynamically shaped streamers. In both scenarios, mul- 
tiple planets can achieve low continuum opacities across 
many tens of AU, while also reducing stellar accretion 
more efficiently than single planets. If planets are divert- 
ing accreting material, transition disks such as J 1604- 
2130 may be promising targets for searches for young 
planets; the additional accretion luminosity will improve 
the prospects for detection, as in the case of the planet 
recently imaged within the transition disk of LkCa 15 
flKraus " fc IrelSidl[20TT ). 



5.4. Comparison to other transition disks 

Surveys with the Spitzer Space Telescope have shown 
that ~ 10 — 20% of circumstellar disks in nearby star 
forming regions have mid-infrared SED decrements in- 
dicative of an inner hole. Most of these objects have low 
millimeter fluxes and are either debris disks or in the late 
stages of uv-photoevaporation ( |Cieza et al.|2008| ). About 
20% have properties that are consistent with ongoing gi- 
ant planet formation (Cieza et al. 2012), as we suspect 
to be the case with J1604-2130. Relative to the full pro- 
toplanetary disk population, therefore, such objects are 
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rare. 



High resolution millimeter interferometry of these po- 
tential planet forming transition disks has resolved the 
inner holes in about a dozen cases to date. Their masses 
range from 0.006 — 0.12 M and inner radii vary between 
15 — 73 AU with no appa rent correlation ( | Andrews et al.| 
2011||Brown et al.||2009| ). J1604-2130 lies below the me- 



dian in terms oi mass but has one of the largest cavities 
yet found. Further, it is one of a handful that are known 
to exhibit strong line emission in the outer disk (e.g., 
|Hughes et al.|2009| ). In a future paper, we will present the 
rich chemistry that we have found here, most likely due 
to the release of previously deplet ed species from grain 
surfaces at the exposed inner rim (jCleeves et al. 2011). 
It also stands out as having negligible stellar accretion, 
which may indicate massive and/or multiple giant plan- 
ets, and which may also maintain the large outer disk 
reservoir over the relatively advanced, ~ 5Myr, age of 
the system. 



6. SUMMARY 

In this paper, we have presented high resolution SMA 
data of the 880 jam continuum and CO 3-2 line emission 
of the transition disk around J1604-2130, a solar mass 
star in the Upper Scorpius moving group. The images 
show a nearly face-on ring of dust and gas and we have 
determined the dust properties through a least squares 
fit to the infrared-millimeter SED and submillimeter vis- 
ibilities. 

We derive a dust mass of 0.1Mj up . The millimeter 
SED slope is a power law out to the longest observed 
wavelength of 2.6 mm indicating the power law distribu- 
tion of grain sizes extends beyond about 0.8 cm. The 
combined visibility- SED fitting is consistent with these 
values and determines the inner edge of the submillimeter 
dust ring to be 72 AU in radius. The stellar photosphere 
is seen out to about 10 jam in wavelength showing that 
the inner 20 AU around the star is dust free, although 
variability at shorter infrared wavelengths suggests some 
dust may sporadically pass through this region. Mid- 
infrared excesses indicate dust emission from the inter- 
mediate zone 20-72 AU although the radially averaged 
surface density is about an order of magnitude lower than 
the outer regions. 

The presence of millimeter CO emission shows that the 
outer disk is more primordial in nature than debris. The 
gas-to-dust ratio is at least unity and probably substan- 
tially higher. Both the CO image and radially averaged 
visibilities suggest that the gas may extend slightly closer 
to the star than the millimeter emitting dust. 

The relatively high mass of the gas-dust disk, its sharp 
inner edge, and the lack of gas accretion onto the star 
together point to a dynamical origin for the inner hole. 
No companions have been found down to brown dwarf 
masses and we suggest that at least one and perhaps 
several giant planets reside within about 40 AU radius. 

J1604-2130 is the most massive disk remaining in Up- 
per Scorpius and therefore perhaps the last holdover from 
the epoch of giant planet formation. It may well be that 
the particular architecture of its planetary system has 
conspired to trap a substantial reservoir of gas and dust 
at large radii, thereby enabling this case study of the late 
stages of primordial disk evolution. 
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